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Abstract. - We study the survival probabihty and the corresponding first passage time density 
of fractional Brownian motion confined to a two-dimensional open wedge domain with absorbing 
boundaries. By analytical arguments and numerical simulation we show that in the long time limit 
the first passage time density scales as pQ{t) ~ i-i+'^(2if-2)/(2e) terms of the Hurst exponent 
H and the wedge angle G. We discuss this scaling behaviour in connection with the reaction 
kinetics of FBM particles in a one-dimensional domain. 



Introduction. — The first passage defines the mo- 
ment at which a random quantity crosses a given thresh- 
old value for the first time. For instance, this could be the 
time when a random walker first leaves a defined region 
of space. Other examples include phenomena as diverse 
as diffusion-limited reactions [Tj, cyclisation of a polymer 
stock market dynamics [3], or search problems [S]. 
The concept of first passage statistics is of quite ubiquitous 
importance in the theory of stochastic processes [6] . 

First passage characteristics have been investigated for 
various geometries. Thus, on semi-infinite domains the 
probability density of the first passage generally displays 
a power-law tail, and the associated mean first passage 
time diverges |^. Recently, remarkable progress has been 
achieved, relating the underlying geometry to the resulting 
first passage behaviour [7]. 

In the case of normal diffusion, characterised by the 
linear growth (r(t)^) ~ t of the ensemble-averaged mean 
squared displacement, the associated first passage prob- 
lems are obtained by solving the associated diffusion equa- 
tion with given boundary conditions or, for simpler geome- 
tries, by the method of images [6]. For anomalous diffusion 
with the non-linear growth [3] 



of the mean squared displacement, no general theory ex- 
ists to obtain the corresponding first passage time statis- 
tics. Only for certain classes of anomalous diffusion mod- 
els analytical approaches for the first passage properties 
are known. Examples include the continuous time ran- 
dom walk model |10| governed by a heavy-tailed wait- 
ing time distribution with diverging characteristic waiting 
time [TTmi] or by a power-law jump length distribution 
with diverging variance (Levy flights) [15]; as well as the 
diffusion on fractal media [5]. Importantly the first pas- 
sage time behaviour differs between different stochastic 
processes despite sharing the same form ([T|) of the mean 
squared displacement with dynamic exponent 211. 

Fractional Brownian motion (FBM), originally intro- 
duced by Kolmogorov [16] and later re-discovered by Man- 
delbrot and van Ness [17], is a generalised Gaussian pro- 
cess, whose increments in one spatial dimension. 



dx{t) ^ i"{t)dt, 



(2) 



arc stationary and normally distributed but dependent. 
Here the quantity [t) represents fractional Gaussian 
noise with zero mean {{^^ (t)) = 0) and autocorrelation 
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Kh of dimension [A'^] = cm^/scc^^ is the anomalous 
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Fig. 1: (a) Stochastic motion of a particle confined to a two- 
dimensional wedge of opening angle O. The particle motion 
motion at time t is described by r{t) = {r{t),8{t)) in polar 
coordinates. The excursion is terminated when the particle hits 
one of the absorbing walls located at = or 0, the moment 
of first passage, (b) Diffusion-limited reaction of three particles 
that move according to FBM in one-dimensional space. 



diffusion coefficient. The fractional Gaussian noise is neg- 
atively (antipersistence) or positively (persistence) corre- 
lated for the cases of subdiffusion and superdiffusion, re- 
spectively. FBM has been widely used to describe anoma- 
lous diffusion phenomena including annual river discharge 
[IH] , stock market dynamics [12] , the motion of tracer par- 
ticles in crowded environments [20^55] . or single-file dif- 
fusion [23]. Despite its popularity and deceivingly sim- 
ple definition in terms of fractional Gaussian noise the 
exact stochastic properties of FBM are not well under- 
stood. The first passage properties are known only for 
the one-dimensional semi-infinite domain. The associated 
long-time asymptotic form of the first passage time den- 
sity p{t) ~ t^~^ was conjectured by Ding and Yang from 
scaling arguments |25j , and rigorously proved by Molchan 
[55]. Recently, it was shown that the Wilemskii-Fixmann 
approximation, originally introduced for the description 
of polymer cyclisation, produces the first passage time be- 
haviour and barrier crossings of particles driven by frac- 
tional Gaussian noise [27]. The difficulty to analytically 
access FBM's first passage properties is related to the fact 
that FBM is a strongly non-Markovian process, which does 
not fall into the class of semimartingales [55] . 

Here we analyse the first passage time behaviour of two- 
dimensional FBM confined to a wedge domain of opening 
angle Q (Fig. [T]). This problem is of particular interest 
as the reaction kinetics of three diffusive particles in one- 
dimensional space can be mapped on this problem [SJUn] ■ 
In the case of ordinary Brownian motion {H = 1/2), the 
wedge problem can be solved exactly using the Green's 
function formalism [5] . Based on analytical and numerical 
arguments we conjecture the first passage behaviour of 
FBM for arbitrary Hurst exponent in a wedge geometry. 



In the following Section we define multi-dimensional 
FBM and the first passage in a wedge domain. We then re- 
view the Green's function formalism for the Brownian case 
and obtain analytical forms for the first passage time quan- 
tities. We proceed to present our main results, the asymp- 
totic scaling forms of the first passage time behaviour as 
function of H and Q. These are numerically confirmed 
in the subsequent Section. We conclude with a discussion 
related to reaction kinetics of FBM particles. 

Fractional Brownian motion in wedge domain. 

Consider the random walk of a particle confined to a two- 
dimensional wedge domain of opening angle 0, see Fig. [TJ 
The particle starts at a position rg inside the wedge and 
moves within the domain until eventually it hits one of the 
absorbing boundaries at a time t for the first time. For 
this process, we obtain the survival probability 



^e(i) 



G{r,t)dr 



(4) 



w 



on the wedge domain W, and the first passage time density 



peit) = -- 



dt 



(5) 



In Eq. (|3]), G{r,t) is the probability density to find the 
particle at position r at time t. Due to the presence of 
the absorbing boundaries, G(r, t) is not normalised, as 
expressed by the survival probability, which decays from 
J5^e(0) = 1 to limt^ao y&{t) - 0. 

Two-dimensional FBM is defined as a superposition of 
independent FBMs for each Cartesian coordinate [30H32] : 
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(6) 



where Xi is the unit vector in Cartesian direction i [i = 
x,y), and ^/^(i) is the fractional Gaussian noise. Due to 
this definition it is clear that for a quadrant geometry 
with 9 = 7r/2 the absorption to either wall decouples. 
We now first address the case of normal Brownian motion 
{H = 1/2), for which exact results for the survival and 
first passage distributions can be found. 

First passage process for H = 1/2. — For normal 
Brownian motion, the diffusion in a wedge domain is de- 
scribed in terms of the Green's function G{r,9;t), that 
satisfies the diffusion equation in polar coordinates, 

d / 1 d 1 \ 



dr' 



r dr d(fi 



(7) 

The associated boundary conditions are G(r, 0; t) ~ 
G(r,Q;t) ~ 0, representing the absorbing walls. The so- 
lution G is completely specified by the initial condition 
G{r,0;O) = 5{r - ro)S{e - 9o)/ro. At long times i, the 
solution G can be approximated by [5] 
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where /^(z) is the modified Bessel function of the first 
kind, which can be expressed by the series expansion 



(zV4)^ 



k=0 



(9) 



The survival probability @ in polar coordinates is then 

re /-oo 



Jo 



rG{r,9;t)drd0, 



(10) 



from which the property ^e(O) = 1 follows directly from 
the sharp initial condition. Using the approximations 
/^(z) « (z/2)'^/r(l + ly) and e-''o/*'^i/2* ~ i for the long- 
time limit in Eq. (j8]), one obtains the scaling expressions 



^e(i) 
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Notably, the wedge angle 6 enters the scaling exponents 
inverse-proportionally; survival and first passage distribu- 
tions decay faster for decreasing wedge angle, as it should. 
We observe an interesting crossover as function of 0: as 
long as Q > it/ 2 the mean first passage time 



r = 



tpe{t)dt 



(12) 



diverges, T — !> oo. In particular, in the half-space limit 
G = TT, the process reduces to the first passage in a one- 
dimensional, semi- infinite geometry (along the y-axis), 
and we find the usual Levy-Smirnov scaling p(t) ~ 
When the wedge angle is smaller, 8 < 7r/2, the mean first 
passage time is finite, reflecting the much higher proba- 
bility to hit one of the two walls of the wedge. Thus, for 
8 ~ 7r/2, the absorption to either of the two walls can 
be viewed as two decoupled Brownian random walks in x 
and y direction. The resulting survival probability is then 
given as the product of two one-dimensional survival func- 
tions corresponding to 8 = tt: o5^^/2(i) = ^7r(i)^ — t~^, 
for which the mean first passage time is logarithmically 
divergent. Remarkably, the scaling relation 



holds more generally for any 8, see relations ([TT|) . 



(13) 



First passage process for H ^ 1/2. — Let us now 
address the case of general Hurst exponent, H ^ 1/2, for 
which no analogue to the above Green's function method 
is known. While there exist dynamic equations for FBM 
in literature BSIIMI of the form 



(14) 



with the time-dependent diffusion coefficient K{t) 



i^^^^, this description cannot fully specify the behaviour 
of FBM in the presence of non-natural boundary condi- 
tions. To see this, we follow the procedure of the preceding 
Section, and find the first passage behaviours 



p%{t) oc t-i-^/e. 



(15a) 
(15b) 



While, naturally, the Brownian case i7 = 1/2 is consistent 
with the results of the previous Section, the results ((T5|) are 
inconsistent with Molchan's result p{t) ^ t^~^ for a com- 
pletely open wedge, 8 = tt. Therefore, for arbitrary wedge 
angle 8 and Hurst exponent H the distributions ([T5|) can- 
not be correct. In fact, the same inconsistency is found 
when one naively applies the method of images to the free 
space solution {ATrKnt^")'^^^ cxp (-xV [AKnt^"])- 

We now argue in favour of a conjecture for the correct 
scaling forms of the first passage time quantities. Let us 
start by recollecting the known special cases, for which 
analytical results are available: 

(i) When 8 = tt, the first passage time quantities must 
converge to Molchan's result for a one-dimensional semi- 
infinite domain 1261: 



P7r(i) 



H-1 



H-2 



(16a) 
(16b) 



valid in the long time limit t ^ 1. 

(ii) When 8 = 7r/2, due to the independence of x and 
y motion the survival probability, by above argument, is 
given as the product of two one-dimensional survival prob- 
abilities, 



it/2 



and thus 



P-R/2it) 



2H-2 



2H-3 



(17a) 

(17b) 



(iii) For H = 1/2, we know the analytical results ([TT|) 
with the full 8-dependence. 

Assuming that the first passage time exponents are sim- 
ple combinations of the Hurst exponent H and the wedge 
angle 8, to satisfy the above three special cases we con- 
jecture the unique scaling form 



Peit) 



t 



7r(2H-2)/(2e) 
-l+7r(2_H"-2)/(2e) 



(18a) 
(18b) 



These results imply that the survival decays faster for 
slower diffusion (i.e., lower value of i?), in analogy with 
Molchan's result and findings for the barrier crossing of 
FBM [27] . Note that the inconsistent results (|15p based on 
the dynamic equation (|14l) follow from our results through 
the substitution iJ — !> 1 — iJ. 

Simulations. — To numerically verify the conjecture 
([T5| for the survival probability and the first passage time 
density, we performed simulations for various wedge an- 
gles and Hurst exponents, and compared the results with 
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Fig. 2: Survival probability =5^6) (i) for wedge angles = tt 
(Left) and Q = n/2 (Right) with the Hurst exponents H = 1/4, 
1/2, and 3/4 (from bottom to top). The lines represent the 
expected scaling behaviours, Eqs. H16ap and (|17ap . 

Eqs. p^ . For this comparison, we here focus on the sur- 
vival probability, which generally shows better statistics. 

Our simulations procedure is as follows. A wedge of 
opening angle Q is constructed by imposing two absorb- 
ing walls at y = and y = (tan &)x in the x-y plane 
(Fig. [Ija)). At i = 0, the particle is located at ro = 
(ro cos ^, ro sin ^) on the line bisecting the wedge. For 
t > 0, the particle undergoes two-dimensional FBM in 
discrete time indexed by n, following the rule 

r(i„) = 

^^'(*™) + '^0 '^"^ ^s^'l^™) + ro sin - .(19) 

\m— 1 ni—1 / 

Here, the two independent fractional Gaussian noise func- 
tions £,i^{t) and C^(i) are generated by the Hosking 
method [3TJ[3S]. The simulation is terminated when the 
particle escapes the wedge domain for the first time. This 
is the survival time for the particle in the given run. From 
repeated runs we construct the survival probability S^q (t) 
as the accumulated number of particles surviving until 
time t. For each given wedge angle and Hurst exponent, 
the survival probability ^e(0 obtained from 100,000 
runs and for initial radius ro = 0.25. 

To confirm that our simulations procedure produces cor- 
rect results for the survival probability, we first consider 
the two special cases of wedge angles Q — n and tt/2, and 
compare them to the predicted scaling behaviours (jl6ap 
and (|17ap . Fig. [2] shows on the left the distributions of 
S^eit) obtained from simulations of a fully open wedge, 
8 = TT, for Hurst exponents H = 1/4 (subdifFusion), 
1/2 (normal diffusion), and 3/4 (superdiffusion). Clearly, 
the distributions follow the predicted scaling behaviours 
for this case corresponding to the one-dimensional, semi- 
infinite domain. Similarly, for the case of a rectangular 



wedge (9 = 7r/2) on the right of Fig. [2] we find excellent 
verification of our simulations method, compared to the 
predicted behaviour o5^7r/2(0 ^ t^^^^ for all values of H. 

We now investigate whether our conjecture (|18ap for 
the survival probability also holds for other wedge angles. 
In our simulation, the survival probabilities for wedge an- 
gles 8 = 7r/4, 37r/4, 57r/4, and Ttt/A were investigated, 
for Hurst exponents H = 1/4, 1/2, and 3/4. As demon- 
strated in Fig. [3l the simulated survival probabilities re- 
markably well follow our conjecture (|18ap for all wedge 
angles and Hurst exponents. Particularly, the scaling ex- 
ponents of the survival indeed decreases for growing values 
of the Hurst exponent H, in line with Molchan's result and 
findings for the barrier crossing behaviour of FBM. 

Discussion. — We studied the first passage proper- 
ties of two-dimensional FBM confined to a wedge domain 
with absorbing boundaries. Starting from special cases, 
for which exact long-time scaling expressions can be ana- 
lytically derived, we obtained functional dependencies of 
first passage time exponents on the Hurst exponent H and 
wedge angle 8. By numerical analysis we confirmed this 
conjecture. We believe that our results are an important 
additional brick in the construction of a clear picture of the 
elusive FBM process. The obtained forms ([T8| of the first 
passage time quantities imply that the mean first passage 
time of a particle from a wedge domain across the domain 
walls, diverges for wedge angles 8 > 8c. The critical 
angle corresponds to 8c = 7r(l — H). 

From the above results one may obtain some insight into 
the reaction kinetics of three particles in one-dimensional 
space, driven by FBM (Fig. IHb)): (1) Surrounded prey. 
Let us consider the diffusion-limited reaction of the cen- 
tral particle with either of the two comrade particles. The 
problem can be transformed such that one particle is con- 
fined between two moving absorbing boundaries [6]. If 
Xi(t) denotes the position of the ith particle at time t, the 
survival condition of the central particle is xi(t) < X2{t) 
and X2{t) < X3{t) for particles sharing the same diffusion 
properties {H and Kh). We regard Xi{t) as the ith. coordi- 
nate of a single particle, that diffuses in three-dimensional 
space. Then, the survival criterion is understood such that 
this particle diffuses inside a domain limited by absorbing 
planes at .ti = X2 and X3 = X2- This is equivalent to 
the first passage problem for two-dimensional motion in 
a wedge domain, the wedge angle given by the two in- 
tersecting planes, i.e., 8mid = 7r/3 (see Ref. [6]). Thus, 
from the distribution ()18ap we find that the survival prob- 
ability for the central particle has the asymptotic form 
=5^mid(i) t^''"-^\ Contrary to our naive expectation, 
the above scaling exponent shows that superdiffusive par- 
ticles (with H > 1/2) survive longer than ordinary Brow- 
nian particles (H = 1/2), and subdifFusive particles have 
a better chance to meet each other. These characteris- 
tics will lead to significant differences in the total reaction 
amount, compare Figs. [2] and [3] Intriguingly the mean 
reaction time ^ is finite for H < H* =2/3, meaning that 
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Fig. 3: Survival probability =5^e(i) for wedge angles Q = 7r/4, 37r/4, 57r/4, and 77r/4, with Hurst exponents H — 1/4, 1/2, and 
3/4 (from bottom to top). The corresponding expected scaling behaviours (cf. Eq. (|18a[) ) are depicted by the full lines. 



reactions always occur for subdiffusive motion, similar to 
the considerations in Ref. [55]. (2) Chased prey. For the 
case that one of the corner particles {xs , say) is chased by 
the two other diffusing particles (xi and X2), the problem 
is reformulated to the task of finding the survival probabil- 
ity for the wedge domain of angle 6ond = 27r/3. Hence for 
identical particles the probability that one corner particle 
survives until time t scales as J^end{t) ~ t3(-H'-i)/2^ ^j^jg 
case the reaction is slower than that of the surrounded par- 
ticle. Note the scaling relation S^mid{t) — =^end(^)' such 
that the reaction of a particle surrounded by two others 
corresponds to the product of two independent corner par- 
ticle reactions. 

We note that closer inspection of the simulated sur- 
vival probabilities shows that in the non-Markovian case 
{H ^ 1/2) the first passage process exhibits a somewhat 
intricate scaling behaviour: the long-time behaviours (jl8|) 
are preceded by another distinct scaling law at intermedi- 
ate times. This behaviour is particularly noticeable for the 
cases of wedge angle < tt, for which first passage events 
occur easily at short and intermediate times. The appar- 
ent intermediate scaling behaves differently, depending on 
the diffusion characteristics. For subdiffusion, the inter- 
mediate scaling exponent is smaller than the long-time 
exponent §(1 — H) and larger for superdiffusion. 

The first passage behaviour displayed by FBM is differ- 
ent from results derived from the diffusion equation with 
time-dependent diffusion coefficient. Latter would lead to 
an inconsistent _ff-dependence, corresponding to the re- 
placement H ^ 1 — H in the correct forms (fT8|) . Similar 
inconsistency occurs if the method of images were applied 
to construct the solution in the presence of non-natural 
boundary conditions. We note in passing that the sub- 
stitution if — > 1 — H m our results for = tt leads to 
the scaling for the first passage behaviour found in the 



analysis of a generalised Langevin equation in Ref. |37| . 

The somewhat counterintuitive behaviour that smaller 
H implies faster decay of the first passage time density, 
may be connected to the fact that FBM is fuelled by ex- 
ternal noise, that is not balanced by friction. FBM, that 
is, does not obey the fluctuation-dissipation relation, in 
contrast to the generalised (fractional) Langevin equation 
[33l[34l[38] that describes FBM in conditions close to ther- 
mal equilibrium, when the fluctuation-dissipation theorem 
is valid. We will investigate the latter behaviour in a sep- 
arate work. 

We flnally note that while for a compact process in one 
dimension the first passage is identical to the first arrival 
to a given position, these processes are no longer necessar- 
ily equal to each other in higher dimensions. This may also 
be the reason why the Wilemskii-Fixman approximation 
reproduces Molchan's result p{t) ~ t^~^ in one dimen- 
sion, but delivers a different result in higher dimensions. 

* * * 

We thank Michael Lomholt and Igor Sokolov for help- 
ful discussion. Financial support from the Academy of 
Finland (FiDiPro scheme) and the European Commission 
through MC IIF Grant No. 219966 LeFrac is gratefully 
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Appendix. — The survival probabilities presented in 
Fig. [3] were obtained with the same initial starting point 
and number of simulation runs for consistency. Depend- 
ing on the wedge angle and Hurst exponent, however, some 
cases are more problematic to show reliable long-time scal- 
ing properties within the time window we used. Here we 
present supplementary results for these cases. 

(i) Fig.|4]depicts the case for 9 = 7r/4 with ro = 2.0 (i.e., 
eight times the value used in the main text) and 250,000 
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Fig. 4: (Left) Survival probability ^e(i) for wedge angle — 
7r/4 and Hurst exponents H = 1/4 (below) and 1/2 (above). 
Here ro = 2.0, and the number of simulation runs is 250,000. 
(Right) Survival probability ,i^e(i) at wedge angles O — Sir /A 
for the Hurst exponent H = 1/4. Here ro = 0.10, and number 
of simulation runs is 100, 000. The solid lines are the expected 
power-laws predicted by Eq. (|18a|) . 

runs. Now the survival probabilities for H = 1/4 and 1/2 
exhibit satisfactory long-time scaling. 

(ii) To obtain extended long-time scaling behaviour for 
H = 1/4 and wedge angle Q = Sir/ 4, the initial distance 
from the origin was decreased to ro = 0.10. The obtained 
result in Fig. 2] shows good agreement with Eq. (|18ap . 
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